MODELING MULTIPLE RISKS: HIDDEN DOMAIN OF ATTRACTION 
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Abstract. Hidden regular variation is a sub-model of multivariate regular variation and facilitates 
accurate estimation of joint tail probabilities. We generalize the model of hidden regular variation 
to what we call hidden domain of attraction. We exhibit examples that illustrate the need for a 
more general model and discuss detection and estimation techniques. 

1. Introduction 

Tail probabilities, especially joint tail probabilities, provide useful risk measures for many ap- 
plications including finance (l3l | , environmental protection [13] and hydrology 0] ■ Multivariate 
extreme value theory (MEVT) is a tool to approximate such tail probabilities but in many common 
circumstances the tool gives an incorrect tail probability approximation of 0. This paper points 
out that even when hidden regular variation (HRV) is not applicable, a more general concept called 
hidden domain of attraction may yield a fix. 

The joint distribution H{-) of a bivariate random vector X = belongs to the maximal 

domain of attraction of a bivariate distribution G{-) if there exist scaling and centering constants 

> and 6^, i = 1, 2, such that for all continuity points x = {x^,x'^) of G, 

(1.1) hm [H{aix^ + blalx' + bl)]'' = Gix\x') 

and both the marginal distributions of G(-), G^(-) and G'^(-), are non-degenerate extreme value 
distributions @, page 208]. The convergence relation (jl.lh is equivalent to the condition that as 
n oo, 

'X^-bl. X^-bl' 



(1.2) nP 



in M+(E), where E = J-oo, oo]^ \ {(-oo, -oo)} or E = JO, oo]^ \ {(0, 0)} or E = J-oo, oo] x [0, oo] \ 
{(— oo, 0)} or E = [0, oo] x J-oo, oo] \ {(0, — oo)} depending on the case. Also, M+(E) denotes the 
set of Radon measures on E and A denotes vague convergence. The limit measure in ()1.2p is 
related to the limit distribution G(-) in (jl.ip as follows: for x = (x^,2;^) E E, 

(1.3) z/({(z\ z^) G E : < x\ < x^Y) = - log(G(x\ x^)). 

Assuming [X^ > u, X"^ > -y] is a rare event, that is, that u and v are sufficiently large, we use 
MEVT to approximate the joint tail probability P{X^ > u, X"^ > v) as 



(1.4) P{X^ >u,X'^ >v)^ -u 



1 f f u — bl, 



n \\ a}. 



, oo 



6! 



.2 



, oo 



However, in the presence of asymptotic independence 0, page 226], (|1.4p approximates the joint 
tail probability P{X^ > n, > v) as zero. Perhaps this approximation is crude and a better 
estimate is possible. 



Key words and phrases. Regular variation, maximal domain of attraction, spectral measure, risk sets. 
S. Resnick and A. Mitra were partially supported by ARO Contract W911NF-10- 1-0289 at Cornell University. 

1 



2 



A. MITRA AND S. I. RESNICK 



If ([O]) holds with E = [0, oo]2\{(0, 0)}, = bl = and some a^, t OO) we obtain multivariate 
regular variation (MRV). If and X'^ are also asymptotically independent, we may improve the 
app roximation of joint tail probabilities if hidden regular variation (HRV) is p resent; see Resnick 



jl4l |. Mitra and Resnick [1^, Resnick jl5l | and the seminal Ledford and Tawn jld.llll|. However, HRV 
requires the distribution of X^ A to have a regularly varying tail and this may not be the case. 
Perhaps X^ A X"^ has a distribution in some maximal domain of attraction other than the heavy 
tailed domain. In this case, HRV cannot be applied to improve joint tail probability approximation 
but the deficiency can be remedied by a more general approach which we call hidden domain of 
attraction (HDA). HRV is a special case of HDA. 

If the distribution of X does not have MRV but (jl.2p still holds, we may retrieve the MRV setup 
by transforming the components of X to C/^(X^)) , where [/*(•) = 1/(1 — and H^{-) 

is the distribution of X*, i = 1,2 161, page 265]. If X^ and X^ are asymptotically independent, so 
are (C/^(X^), C/^(X^)) and assuming U^{X^) A C/^(X^) has a regularly varying tail, we may seek 
HRV. Statistically this is problematic since we do not know [/*(•), i = 1, 2. This can be dealt with 
in various ways, none of which is completely satisfying or easy and a potential advantage of the 
notion of HDA is that it does not require that we transform components. 



1.1. Outline. Section [1 .21 reviews frequently used notation. In Section [2l we define hidden domain 
of attraction for the standard case, when both the components of the risk vector have the same 
distribution. Section [3] deals with HDA in the non-standard case, where we drop the identical 
distribution assumption for X^,X^. In both Sections [2] and [3l we exhibit examples which satisfy 
our model and discuss estimation procedures of limit measures that appear in the limit relations 
of the model. Section U] discusses the detection techniques for HDA and estimation of joint and 
marginal tail probabilities. We conclude with a few remarks in Section [5j 



1.2. Notation. For simplicity, this paper is restricted to two dimensions. For denoting a vector 
and its components, we use: 



X = [x ,x ), X* = i-th component of a;, i = 1, 2. 

Multivariate intervals or rectangles are denoted {x,y],[x,y], etc where, for instance, {x,y] = 
{x^,y^] X The vectors of all zeros, all ones and all infinities are denoted by = (0,0), 

1 = (1,1) and oo = (00,00) respectively. We write x*-^-* = V x^, x*-^^ = x^ A x^. So, the 
superscripts denote components of a vector and the ordered component is denoted by a parenthesis 
in the superscript. 

We express vague convergence [l^, page 173] of Radon measures as A and weak convergence of 
probability measures 0, page 14] as Denote a point measure with points {xj} in a nice space 
F by where for x G F and B CF, 



1, if X G B, 
0, ifxG^^ 



Write M_|_(F) for the set of non- negative Radon measures on a space F topologized by the vague 
topology. 

For a one dimensional distribution -F(x), set F := 1 — F. The inverse of a non-decreasing function 
ip{x) is ijj^{x). 
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2. Standard case hidden domain of attraction 

Suppose that a bivariate random vector X = with distribution H{x) belongs to the 

maximal domain of attraction of an extreme value distribution G 16, page 265], = X"^ and X^ 
and X"^ are asymptotically independent so that (jl.ip is satisfied with 6^ = 6^ = 6„, = = a„> 
and G{x,y) = G^{x)G^{y) for x,y e M, n G N. Thus, (fO]) becomes 



(2.1) nP 



X-bnl 



an 

in M+(E), where E = [-oo, oo]\{-oo} or E = [0,oo]\{0}. Since G{x,y) = G^{x)G^{y) for 
j;, y € M, the relation of and G given in (II. 3p gives for £C € E, 

v{{z eE: <x\z'^ < x'^}^) =-\ogG^{x^) + -logG^{x^) 

(2.2) =i/({z eE: < x^}^) + i/({z G E : < ^2|C^^ 

The standard case contains the additional assumption that X^ = X"^, which reduces (jl.2p to (j2.ip 
and reduces possible choices for E. The cone E = [0, oo]\{0} is chosen only when H has MRV. 
From ()2.ip . the maximal component of X satisfies as n — > oo, 

(2.3) nP(X(i) > any + fe„) ^ i^({z G E : z^^) > y}), {y, y) G E, 

so X^^^ is in a maximal domain of attraction of an extreme value distribution and the distribution 
of X*-^-* characterizes a„ and 6„ given in (j2.ip . Using one-dimensional extreme value theory, we can 
choose a„ and 6„ in such a way that 

(2.4) i^iy) := u{{z G E : > y}) 
takes one of the following forms: 

^(^) = 1 0, otherwis;, ^ > 

^p{y) = e^, if 7 = 0, 

^(^) = | (-y)-V7, otherwise, ^ < °' 

where 7 is the extreme value index of the distribution of [H, page 9]. This remams our 

standing assumption for the following discussion. 

We define a sub-model of MEVT called (standard case) hidden domain of attraction (HDA) . HDA 
helps approximate joint tail probabilities in the presence of asymptotic independence and includes 
HRV as a special case. If E is either [— 00, cx3]\{— 00} or [0,oo]\{0}, define E° as (— oo,cx)] or 
(0, 00] but see Remark 12.2112]) before jumping to erroneous conclusions that E'^ C E is always true. 

Definition 2.1. The distribution of X = (X^, X'^) has standard case hidden domain of attraction 

on the cone E° if (i) X^ = X^\ (ii) (12. ip and (12. 2p hold; (iii) there exist positive scaling and real 
centering constants {cn} and {dn} and a non-zero measure G M+(E'') such that in M+(E'^), 

(2.5) nP[(X -d„l)/c„ G •] A i/°(-) (n^oo). 

We emphasize that the definition requires that X^ = X"^ and that the distribution of X belongs 
to the maximal domain of attraction of an extreme value product measure G with exponent measure 
V. Some other remarks: 
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Remark 2.2. (1) Hidden regular variation assumes (j2.ip is satisfied on the cone E = [0, oo]\{0} 
with bn = and a„ t c« and (j2.5p is satisfied on the cone E*^ = (0, oo] with dn = and 
c„ t oo. Moreover, an/cn — )• oo as n ^ oo. Hidden regular variation is a special case of 
hidden domain of attraction. HRV is the only sub-model of HDA where the cone E*^ in 
(f23D is E° = (0,cx)]. 

(2) From (j2.5p the minimum component of X satisfies, 

(2.6) nP[X(2)>c„y + d„]^z.0((y,oo]x(y,oo]) {{y,y) eE^), 

and therefore, the distribution of X^^^ belongs to the maximal domain of attraction of an 
extreme value distribution 0, page 4]. When HRV exists, the distribution of X^^^ has a 
regularly varying tail and is hence in the domain of attraction of the Prechet distribution. 
HDA allows the additional cases where the distribution of X^'^^ belongs to the domain of 
attraction of the Gumbel or the Weibull distribution. 

The distribution of X^"^^ determines the scaling and centering constants {cn} and {dn} 
and the cone E°. As illustrated by Example O even if E = [0, oo] \ {0}, the cone E° could 
be {—oo, oo] and E'^ is not necessarily a sub-cone of E, as was the case for HRV (l^ . 

(3) From (|2.3p . we get that X^^^ belongs to the maximal domain of attraction of some extreme 
value distribution. Since X« >X(2), the convergence relation ()2.ip puts some restriction 
as to what possible convergences can hold in ()2.5p . For example, if (12. ip is satisfied with 
X^^^ being in the Gumbel domain of attraction, then HRV can never hold on E'^ since the 
tail of X^^) cannot be heavier than the tail of X^^\ 

2.1. Semi-parametric structure of u^. The limit measure in (|2.5p has a semi-parametric 



structure which assists estimation (as in [121 ] for HRV) and which characterizes the class of possible 
limit measures as a class indexed by a real parameter and a set of probability measures. 

To understand this semi-parametric structure, proceed as follows. Let H^'^\-) be the distribution 
of X(^) and define the function as 

(2.7) ^^y):=[z.O((y,oo]x(y,oo])]-\ 
where I'^i-) is given in ()2.5p . From ()2.6p we get 



(2.8) nH(^){cny + dn)^[r{y)] (y € M). 

Hence, from the one-dimensional extreme value theory, H^'^'^ is in a maximal domain of attraction. 
Let 7*^ be the extreme value index of if . Assuming Cn and dn are chosen suitably [igI, page 9] , 
tp^i') must take one of the following three forms: 

(2.9) V'°(y) = eJ', if7° = 0, 

ip^iy) = 1 2/ > 0, if ^0 ^ Q 

\ (— y)"^^'^ , otherwise. 



Henceforth assume that c„ and dn are chosen so that (|2.9p is true. Define 

(2.10) C/(2)(x) = 1/(1 -F(2)(x)), 

and the following helps us identify the semi-parametric structure of i^^. 
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Proposition 2.3. The convergence in ()2.5p that defines HDA is equivalent to the regular variation 
on (0, ool, 



(2.11) 



nP 



G 



Az>°(.) 



(C/(2)(X1),C/(2)(X2) 

where U^'^\-) is defined in (|2.1U|) and is a Radon measure on (0, oo] that is related to the limit 
measure in p.Sp by 

(2.12) i>^{{x\oo] X {x\oo]) = i/0(((V°)^(x^),oo] x ((V'°)^(x2), oo] ) . 

The measure satisfies the scaling property: 

(2.13) pO(c.) = c-^P°(-), oO. 



(a; G (0, oo]). 



Remark 2.4. (i) Proposition 12.51 below shows that the hmit measure u'^ is determined by a 
probabiUty measure S^. Thus the family of limits in (I2.11|) is indexed by probability measures 
and Proposition 12.31 shows that has semi-parametric structure: the probability measure 
determines and given 7^, we get V°(-) from ([22]) and then applying ([232]) . we get . 
(ii) If the support of the distribution of X^'^^ is smaller than that of X*, i = 1,2, then ^/(^^(X*) 
could take the value 00 with positive probability. Hence, in the following discussion, we treat 
i = 1, 2, as extended random variables. 

To see that (|2.1ip imphes ([23]) . observe that for jc G E°, 
X^-dn 



Proof of Proposition \K 
^X^ - dr, 



nP 



> X 



nP 



U^^\X^) Ui^)(cnX^+dn) U^^\X^) Ui^\cnX^ + dn) 



n 



n 



n 



n 



^ P0(«(xl),0o] X (<(x2),CX)]) = l/°((x\cX)] X (x2,00]), 

where the convergence follows from (12. Sp and (12. lip and the last equality follows from (I2.12P and 
the forms of ^'^(•) given in (]2.9p . Hence, (j2.5p holds. The converse is similar and is omitted. □ 

The scaling property (I2.13P implies that we can express (12. lip in an alternate coordinate system 
that transforms the limit measure into a product. From Proposition 12.31 if the distribution of X 
satisfies Definition 12.11 and has HDA, then {U'^'^\X^),U^'^\X'^)) has regular variation on (0,oo]. 
So, using (I22D, and (f232D . we get that 

1 2^ 



pO((l,oc]) 



1. 



This plus the scaling property (j2.13p implies p''([l,oo]) = 1. This, Proposition 3.1 of |12] and 
Proposition 12.31 yield the equivalent convergence in alternate coordinates given in ()2.14p below. 

For Proposition 12.51 we need the following: let vi be a Pareto measure on (0, 00] satisfying 
v^{{y,oo\) = for y > 0. Since [/(2)(.) is non-decreasing, C/(2)(x(2)) = U^^HX^) A U^^\X^). Set 

= {a; G (0,oo]2 : x^^) = 1}. 



Proposition 2.5. 

(2.14) nP 



The convergence in (j2.1ip is equivalent to 

'u^^{x^ f u(^){x') u(^){x^) y 



n 



\^C/(2)(X(2))' C/(2)(X(2)) J 



A Z^i X 50(.), 
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on (0,00] X JN^^), where is a probability measure on The relation between in (j2.1ip 

and is 



(2.15) 

for r > and Borel sets A C 5^^^^) . 

We call the probability measure standardized hidden spectral measure. 

HRV is a special case of HDA. For HRV, a similar spectral measure 5(2) (.) is defined on 
called the hidden spectral measure by 



P 



G . a:(2) > tj ^ 5(2)(.) (t^oo). 



The standardized hidden spectral measure S'^{-) on iJK^^) given in (j2.14p is 



(2.16) 



(C/(2)(X1),?7(2)(X2)) 



> t 



C/(2)(X(2)) 

where the convergence holds as t — )• x^{2) = sup{?/ G M : H^'^\y) < 1}. If HRV exists, Xjj(2) 



00. 



2.2. Examples. We give examples of distributions that possess multivariate regular variation with 
asymptotic independence. Each has hidden domain of attraction but not hidden regular variation 
emphasizing the need for a concept beyond HRV. 

Example 2.6. Suppose, Wi, W2 F{-) and Zi, Z2 D{-), where F and D belong to the maximal 
domains of attraction of the Frechet with index a = 1 and Gumbel distributions respectively. Let 
be a Bernoulli random variable such that P[B = 1] = 0.5 = 1 — P[B = 0]. Assume the random 
variables Wi , W2 , Zi, Z2 and B are mutually independent and define a bivariate random vector X 
as 

X = {X\X^) = B{Wi,Zi) + (1 - B){Z2, W2). 

We show that the distribution of X has MRV. It suffices [3] to verify that tiX^ V t2X'^ has a 
regularly varying tail for any ^1,^2 > 0. Since F has a regulary varying tail and 

P[tiWi V t2Zi > x]r~^P[tiWi > x], (x 00), 

tiX^ V t2X'^ also has a regularly varying tail. Thus for appropriate a„ t 00, 

nP[X/an G •] Az.(-) 

on E = [0,oo]2 \{0}, where u {{{0,x^] x [0,x^]Y) = \ ([x^Y^ + (2;^)"^) • Therefore, the distribu- 
tion of X has MRV with asymptotic independence on E. 

Furthermore, X^^) belongs to the maximal domain of attraction of a Gumbel distribution and 
therefore, HRV does not exist. To see this, without loss of generality [H, assume that D is a 
von-Mises function [lil . page 40] and for specificity assume the right endpoint of D is infinite. The 
form of the tail is 

D[x) = Ce- IMy^r'dy ^'^oo 



where c > is some constant. Likewise, assume without loss of generality 16|, page 58] that 
F satisfies xF' {x) / F{x) 1, where 1 is the index of regular variation of F. Then H^'^\x) := 
F[x)D{x) is the tail of a von-Mises function with auxiliary function 



fm^){t) = tfD{t)/{t + 1 • fD{t)) * ~~ fD{t) 
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Since 



~ t + 



fumit) H(^){t) F{t) ' D{t) ^ ' foit)' 

From (l6l . Corollary 1.7, page 46] i/^^^ belongs to the maximal domain of attraction of the Gumbel 
distribution. 

Next we choose scaling and centering constants {c„} and {dn} in (gS]) so that = e*' (one 

of the forms in (]2.9p ). The usual choices are 0, page 40] dn= [l/H(^)y (n) and c„ = //^(2)((in)- 
Now, observe that for a; € E*^ = (— oo, oo]^, as n ^ oo, 

nP[X^ >dn + CnX^.X^ > dn + CnX^] 



{dn + CnX^) D {dn + Cnx"^) + {dn + C^X^) D {dn + 

F {dn + CnX^) 



n— 



CriX 



n (F {dn + CnX^) \ 
2 [F{dn + CnX^) ) 
1 1 2 



{dn + CnX ) + - [ = 



2 \F{dn + CnX^) 



if(2) {dn + CnX^) 



The convergence follows from the facts that F is regularly varying, Cn/dn and (|2.8p holds with 
ip^{y) = e^. Therefore, as in Definition 12.11 the distribution of X has HDA on E*^ = (—00,00]^ 
with limit measure 1^^ such that for x = (x^,^^) E E'^, 



u'{{x\^]x{x',^]) = -{e 



+ e~ 



Thus the distribution of X is regularly varying on E, has HDA on E*^, but does not have HRV since 
HRV requires the distribution of X^^^ to be in the domain of attraction of the Frechet distribution 
[H, page 54]. 

Example 2.7. Suppose, U ~ Uniform([0, 1]). Define the random vector X as 



X = {x\x^ 



1 



1 



U' l-U 



Now, note that for x^,x'^ > 0, 2n > {x^) + {x'^) 



n(l - P[l/U < 2nx\ 1/(1 - U) < 2nx^]) = n(l - P[U > (2nx^) \ C/ < 1 - (2nx^) ']) 

= n{l - (1 - {2nx^y^ - {2nx^y^)) ^((a^^)"^ + {x^)~^), 

as n — )• 00. Therefore, on E = [0,00]^ \ {0}, nP [X /2n € •] A iy{-) where the limit measure u 
satisfies u {{[0,x^] x [0, x^])"^) = ((x^) + (x^) )/2, for a; € E and thus the distribution of X has 
MRV with asymptotic independence. 

Also note that for {(x^, x^) € ( — 00, oo]^ : x^ + x^ < 0}, and large n. 



„2^-ll 



nP 



(2.17) 



X^ >2 + 



2x 



n 



n + 1 
n + 1 



,X^ >2 + 



2x' 



nP 



n + 1 
n + 1 + 2x2 



n + 1 n + 1 + 2x2 

U < ^—^-r^,U > 2^ + 2 + 2x2 



2n + 2 + 2x1 



2n + 2 + 2x1 2n + 2 + 2x" 



1 

n I — 



X 



1 



X 



2 2n + 2 + 2xi 2 2n + 2 + 2x2 



h{-x') + {-x')), 
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as n — )■ OO. Similar calculations show that for {{x^, x^) G {—oo, oo]^ : + > 0}, 



(2.18) nP 



2x „ 2x 



> 2 + > 2 + 



0, 



n + 1 n + 1^ 

as n — >■ oo. Therefore, the distribution of X has HDA as in Definition O on E° = ( — oo, oo] with 
limit measure such that for a; G E*^, 

i((-xi) + (-x^)), ifx^ + x2<0 
otherwise. 



((x\oo] X (x^,oo]) 



From (j2.17p and (12.181) i t also follows that X^'^^ belongs to the domain of attraction of the reversed 
Weibull distribution page 59]. So, the distribution of X has HDA, but does not have HRV and 
furthermore, is not a subset of E. 

2.3. Estimation. Because the marginal distributions are assumed to be the same, the standard 
case is somewhat unrealistic for applications but it is important to understand estimation for this 
case before moving on to more realistic scenarios. 

To estimate joint tail probabilities, we first estimate the limit measure given in (12. Sp . Let, 
{X, Xi, i = 1,2, ■ ■ ■ , n} be iid with a common distribution satisfying (j2.5p . From (12. Sp . 



1 

(2.19) -ry^,ef Xi-d(n/k)i \{-) t^°(-) (/c ^ oo, ^ oo), 

k ^ V ) 

in M+(EP) [H, page 139]. From (j2.6p . the distribution of X^^^ determines Cn and dn- The iid data 
{X^"^^ : i = 1,2, ■ ■ ■ ,n} allow estimates (0], [H, page 93]) of c{n/k) and d{n/k), denoted by c{n/k) 
and d{n/k), satisfying 

(•2 20) c(n/A:) ^ rf(n/fc) - d{n/k) p_ ^_ 

c{n/k) ' c(n/k) ' 

Therefore, we get the joint convergence 

r9 9n (^S^. d{n/k)-d{n/k) c{n/k) \ p ^ 

in M4.(E'^) X M?. Apply the almost surely continuous map h, a) ^ i^(a[(-) + bl]) in (j2.2ip and 
we get the following proposition: 

Proposition 2.8. Let, {X,Xi,i > 1} he iid with common distribution satisfying (|2.5p . Then, 

■°(-) {k ^ oo, n/k ^ oo). 



1 



m M+(E°). 

Estimation of v^{-) in Proposition 12.81 does not exploit the semi-parametric structure discussed 
in Section [2.11 and has the disadvantages that (a) there is no guarantee the estimator v^i') is even 
a member of the class of possible limit measures; and (b) we are required to estimate c(-) and d{-). 
These disadvantages are overcome using the semi-parametric structure as was done for HRV in 
[l^ . We need to estimate the the extreme value index 7*^ of the distribution of X^"^^ as well as the 
standardized hidden spectral measure Since {xf'^ : i = 1, 2, • • • , n} is iid data, estimating 7'^ 
of is a standard procedure 0, page 65] so we concentrate on estimating S^{-). A modification 
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of a ranks method [1, 0,(3) [3 to obtain an estimator of avoids the need to estimate c(-) and 
d{-). For 2 = 1, 2, • • • , n, define 

(2.23) R]'^'^ = \{j : > X}}\ and Rf^'^ = \{j:Xf>Xf}\, 

where | • | denotes size of a set. Note that < -R^ < n for i = 1, 2, • • • , n, j = 1, 2. Also notice that 
since R]'^^^ V Rf^^^ = \{j : xj^^ > ^}|, 1 < r]'^^^ V i?f < n for i = 1, 2, • • • ,n. Proposition 
12.91 gives an estimator of which we can modify to get an estimator of S^{-). 

Proposition 2.9. We have in M+((0, oo]^), 

1 " 

(2.24) 9^{-):=-Y,eU/jiU^)^,/jiM2)\{-)^i^''{-) {k ^ ^,n/k ^ 

i=l ^ ' 

where (f2lT]l defines and (l2:23D defines r]'^^^ and R^^^K 
Proof. From (j2.5p and definition of ^'^(•) given in (12. 7p . we have in D((0,oo]), 

1 

^E^((xf)-.(nA))/c(n/.))((-'°°]) ^ [^°(-)r'- 
i=l ^ ' 

Hence [H, page 58], inverse functions also converge in distribution in Die/i((0, oo]), the space of 
left continuous functions with finite right limits, 

(2.25) inf{x : ^E^((xf)-.W.))/eK.)) ^ V^} ^ mf{x : [^°(x)]"^ < 1/.} = (^°)-(.). 

Write the order statistics of {X\ , . . . } as X)^^-^ > • • • > ^^^^^ and observe the left side of 
(12:251) is 



n 



(2.26) inf{x : ^ . . /,.^ ((x, oo]) < kjs} 



'X^i^lM)-d{n/k) 



1=1 



(Xr^-d(n/fc))/c(n/fc)j^^-"'-J' - '"J 1 c(n/fc) 



From (f2?T9]) . ([21251) and (fZ^ we get that as A: oo and n/A; oo, 

vfc^ { ^'c{n/'k)°^^ ) ' c{n/k) ' c{n/k) 
(2.27) ^ ^K(-),((V'°)"(«),W"(0)) 

in M4.(E'^) X Dieft{{0, oo]) x DieftdO, oo]). Using the scaling technique as in Resnick [H, page 311] 
we get from (j2.27p that as A: — >■ oo and n/k oo, 

(2.28) 
1 " 

lEl/xW^) x^>x(^) . ^^°(((V'°)"^(s),oo] X ((/)^(t),oo]) =pO((.,oo] X (t,oo]), 
in Ae/t(0, oo]) X Ae/i((0, oo]). Since the left side of (12:281) is 

i=l ^ ' 1=1 ^ ' 

we have proven (|2.24p . □ 
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Proposition I2.9l vields an estimator of the limit measure uix S^{-) and then an estimator of S''^(- 

Proposition 2.10. The convergence in (j2.24p is equivalent to 
^ I " 

(2.29) i,^y,s\{-):=-Y,^f ^ ^^1.(2)^^2,(2) ^1,(2)^^2,(2)^^ (•)^z^ix50(-) 



«=1 I p1>(2)wr2,(2)'1 ^1,(2) . ^2;i2r 



in M+((0,oo] X (Jb^^^^), where vi x S^{-) is given in Proposition W. 

Proof. The proof uses Proposition 12.51 and follows exactly similar steps as that of Proposition 3.7 
of [1^. It is based on the map x ^ {x'^'^\ x / x^"^^) . □ 

From the convergence in (j2.29p . we construct a consistent estimator of S''^(-): 
YJi=i^ f , /„i,(2)^,„2,(2) „i,(2)^„2,(2)\\ ([l,oo] X •) 



V b1'(2)vr2,(2) '1 TJJT) ' 2,(2) 

(2.30) SS(.) := ^' m n ^ 

m M+(5K(2)). Hence, we have obtained a consistent estimator for both the extreme value index 7" 

and the standardized hidden spectral measure S^. 

7 (2) 

It is possible that R- = for some j = 1, 2 and some i = 1, 2, • • • ,n and thus division by zero 
may be indicated in ()2.30p . Though theoretically justified, this is not desirable when writing code 
for an estimator. The continuous bijection T : (5^^^^ 1— )• [0, 1] given by T : a; 1— )• x^/(x^+x^) provides 
an instant remedy. We use the convention that (X3/00 = 1 and l/oo = 0. Using this transformation, 
(|2.30p becomes 

EILi^/ ^ / «i.(2) \\([i,oo]x-) 

I 1,(2) 2,(2) 'I 1,(2)" 2,(2) I I 

(2.31) ^ V" m \^ ^ ^ ° ^ ^) 

fll>(2)^^2,(2) 

in M_|_([0, 1]). Since T is a continuous bijection, we retrieve 5^^ from o T~^. 



3. Non-standard hidden domain of attraction 

To provide more scope for applications, the non-standard case no longer assumes that =X'^. 
However, we have found that to construct a coherent estimation theory requires careful considera- 
tion of the definitions. As in the standard case, the goal is to approximate marginal and joint tail 
probabilities. 

3.1. How to proceed? In order for (|1.2p to hold when X has different marginal distributions, one 
typically needs different centering and scaling constants for the two components of X. Traditional 
theory [l^, page 277, Proposition 5.15] proceeds by standardizing each component. However, 
a theory of hidden domain of attraction that follows this approach encounters problems in the 
estimation procedure that we could not resolve without strong second order conditions. 

We deviate from the traditional MEVT treatment by requiring that both components in (|2.ip 
have the same centering and scaling but permitting the limit measure to have one zero marginal. 
By a zero marginal, we mean that either the limit measure z^(-) in ()2.ip has the property 

i^{{zeE: >y}) =0 ((2/,2/)gE) 
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(3.1) nP 



in M+(E^). 



or the same holds with in place of z^. This could happen, for instance, if the tail of is lighter 
than that of or vice versa. 

In the non-standard case, if we assume ()2.ip . the limit measure i' may satisfy: 

(i) u has a zero second marginal: for {y,y) € E, (|z G E : > y}) = 0; 

(ii) u has a zero first marginal: for (y,y) G E, (j^; G E : > y}) = 0, 

(iii) the cases (i) and (ii) do not hold, but u ({2; £ K : > x, z"^ > y}) = for (x, y) G E, 

(iv) for {x,y) G E, ({z G E : z^ > x, z"^ > y}) > 0. 

Case (iv) means ()2.ip yields non-zero estimates of the marginal and joint tail probabilities, so 
in this case we have no need to define HDA. The definition and analysis of HDA in case (iii) is 
the same as the standard case discussed in Section [2j The definition and analysis of HDA are very 
similar for cases (i) and (ii) so focus only on case (i). 

For our definition of HDA, a relevant state space is E^ where either E^ = [—00, 00] x (—00, 00] 
or E'~' = [0, 00] X (0, 00]. We will see later that for case (i) it is possible to find HDA on both the 
cones E'"' and E*^ in sequence. Thus, compared to the standard case, our estimation procedure here 
might involve analyzing HDA on the additional cone E'~'. 

Definition 3.1. The distribution of X = (X^,X'^) has hidden domain of attraction on the cone 
E^ if dH]) holds with the hmit measure the second marginal of u is a zero measure and in 
addition, there exist constants e„ > and /„ G M and a non-zero measure z^'"' G M+(E^) such that 
as n — >■ 00, 

From (13. Ij) it follows that for {y,y) G E'^, as n — )■ 00, 

(3.2) nP[X^ > eny + /„] ^ z/^ ({(„, v) G E"" : v > y}) . 

Therefore, the distribution of X^, the second component of X, belongs to the maximal domain 
of attraction of an extreme value distribution 0, page 4]. Using one-dimensional extreme value 
theory, il^^{y) := [z^^ ([0, 00] x (y, oo])]~ must take one of the following forms [l6l . page 9]: 

^ [ 0, otherwise, 

(3.3) V"(y) = e^yeM, if 7^ = 0, 

^ (^) = 1 (-y)-V7-, Otherwise, ^ < °- 

The parameter 7'"' in p.3p is the extreme value index of X"^ . We can and always do choose {e„} 
and {/n} in such a way that takes one of the above forms. 

Remark 3.2. We make a few remarks about Definition 13.11 

(1) Since (j2.ip holds, so does (j2.3p and therefore the maximum component X^^^ belongs to the 
maximal domain of attraction of some extreme value distribution. Since > X2, the 
convergence relation ()2.ip constrains the possible convergences in (j3.ip . For example, if 
(|2.3p has X^^^ in the Gumbel domain of attraction, then the distribution of X"^ cannot have 
a regularly varying tail. 

(2) The distribution of X"^ determines the cone E'"' and ()3.2p yields the scaling and centering 
constants {e,„} and {/n}- If the distribution of X"^ is in the Frechet domain of attraction, 
E'~' = [0,00] X (0,00] and otherwise, E'~' = [—00,00] x (—00,00]. 
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There are two possibilities for the hmit measure in (j3.ip : 

(i) the hmit measure puts zero mass on ah sets (x,oo] x (y,oo] for € E^; or 

(ii) the hmit measure puts non-zero mass on one of the sets (x, oo] x (y, co] for G E^. 

The semi-parametric structure of discussed in the next section imphes that for case (ii), 
i^'~'((x,(X)] X (?/,cx3]) > for ah (x, y) G E'^. So, in case (ii), we get non-zero estimates of joint 
tail probabilities and since we accomplished our goal there is no reason to seek further instances 
of HDA. However, in case (i), the measure z^'~'(-) will not provide non-zero estimates of joint tail 
probabilities. A potential solution is that HDA could still exist a smaller cone such as E°. The 
following definition formalizes the concept. For this definition, the state space is E°, where E° = 
(-00,00]^ or E° = (0,oo]2. 

Definition 3.3. The distribution of = (X^,X^) has hidden domain of attraction on the cones 
E^ and E° if Definition O holds, the the hmit measure in (|3.1|) puts zero mass to all sets of the 
form (x,c«] X (y, 00] for (x,y) G E^, and in addition, there exist centering and scaling constants 
{c„} and {dn} and a non-zero measure G M+(E'') such that as n — ?> 00, 



(3.4) nP 



X-dr,.l\ 

G 



i/°(-) in M+(E' 



0\ 



As noted before in (j2.6p . the scaling and centering constants {c„} and {(i„,} in ()3.4p are char- 
acterized by the distribution of X^'^\ the minimum component of X. Recall the definition of 
given in (j2.7p . As was done in the standard case discussion, we choose the scaling and centering 
constants {c„} and {dn} in (|3.4p so that 1/^" takes one of the forms given in (|2.9p . Also, whether E" 
in (j3.4p is (—00,00]^ or (0, cxj]^ is determined by the distribution of X^'^\ 

3.2. Semi-parametric structure of i^^. Both limit measures i^^ of (j3.ip and of (j3.4p have 
semi-parametric structures. Since the semi-parametric structure of i^^ was discussed in Section [2.11 
we concentrate only on the semi-parametric structure of and proceed as follows. 
Recall that the distributions of X and X'^ are H and H^. Define 

(3.5) U\x) = 1/{1- H\x)). 

The following proposition relates p.ip to a regular variation condition on [0, 00] x (0, 00]. Its proof 
is similar to that of Proposition 12.31 and is omitted. 

Proposition 3.4. Convergence in (13. ip is equivalent to regular variation on the cone [0, 00] x (0, 00], 

■c/2(xi) u^x^y 



(3.6) nP 



n n 



z>"(-) (in M+([0,oo] X (0,00]), 



where (|3.5p defines U'^{-) and is a Radon measure on [0, 00] x (0, 00]. The limit measure is 
related to the limit measure in in (|3.ip by the following relation: for (x"*^, x^) G [0, 00] x (0, 00], 

Z>n((x\oo] X (x2,00]) = U''{{{^l,''r{x^),^] X ((^^^^^(^2)^ ^] ^ ^ 

(3.7) ^^([0,xi] X (x2,oo]) = z.n(^{2 eE^ : < {x'), > {i;'')'- {x^)}y 

The measure t''~'(-) satisfies the scaling property: 

(3.8) P^(c-) = c-^i>^(-) 00. 
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Remark 3.5. (i) On the semi-parametric structure of u^: We will see that a probability measure 
S'^ on [0, oo] determines the limit measure i)^. The parameter 7^ and the probability measure 
on [0, 00] determine v^, since given 7'^ and measure S^, we get the function in (j3.3p 

and D'^ which through (j3.7p determines i^'^. 
(ii) If the support of the distribution of is smaller than that of X^, then U'^{X^) could take 
the value 00, so in this case we consider U'^{X^) as an extended random variable. 

The method that shows £''^([1, cx)]^) = 1 also shows t'^([0, 00] x [l,cx)]) = 1. Proposition 4 of 
and Proposition 13.41 give a convergence relation in new coordinates. 



Proposition 3.6. The convergence in (jS.ip is equivalent to 
'C/2(X2) u\X^] 



(3.9) 



nP 



n 



{in M+((0,oo] X [0,00])), 



where ui is a Pareto measure on (0, 00] satisfying z^i((x, 00]) = x ^ for x > 0, and is a probability 
measure on [0,oo], called the standardized hidden spectral measure. The relation between given 
in ([2TT]) and is 

(3.10) {{x G [0, 00] X (0, 00] : > r, x^/x'^ G A}) = r'^S^iA), 

which holds for all r > and all Borel sets A C [0, 00]. 



3.3. Examples. We give examples of distributions oi X = (X^jX^), X^ ^ X"^, and which have 
HDA. In Example 13. 7| the limit measure of ()3.ip puts zero mass on all sets of the form (x^, 00] x 
(x^, 00] for G and HDA also holds on E°. In Example ESI of puts non-zero mass on 
sets of the form (x^, 00] x (x^, 00] for x G E^. 

Example 3.7. Let X^ ~ exp(l), X"^ ~ exp(2) and X^ and X"^ be independent. Then we get 



n 



[l-P [X^ - logn < x\x2 - logn < x^]) 



n 



-(logn+a;-'-) 



-2(log n+x2) 



as n 



00, 



which implies (12. Ij) holds on E = [— 00,00]^ \ {(— cx),oo)} with u (([— (X),x^] x [— oo,x^])'^) = e~^^ 
and u puts mass only on (—00, 00] x {—00} and v has zero second marginal. So we seek HDA on 
E'~'. As n ^ 00, 



nP 



X' 



and also as n 



nP 



X' 



log n 

00, 
log n 



< x^,X^ 



> x\X^ 



logn 2 
„ > X 



logn 2 
— ; — > X 



n 1 



_2(l2|2i+x2) 



ne 



0. 



Thus, HDA exists on E" 



-00, 00 X 



(—00,00] with limit measure u^, where z^'~'([— 00, x^] x 



(2; ,00]) 



and iy^{{x^, 00] X (x^, 00]) = for x G E'~' so z/'"' concentrates on {—00} x (—00, 00]. 



After peeling away both lines through —00, we look for HDA on E*^. A hint for how to proceed is 
provided by X^ A X'^ ~ exp(3). Note that as n — >■ 00, 



nP 



X' 



log n 



> x^,X^ 



logn 2 
— > X 
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Thus HDA exists on E*^ = (—00, 00]^ with limit measure u^, where ^^{{x^, 00] x (x^, 00]) = g-C^^+Sa:^) 
for X e E°. 

For this example, Definition 13.31 holds and HDA holds on both the cones E'"' and E'', 

Example 3.8. Suppose £^1, £'2, -E3 are iid exp(l) random variables independent of i? ~ Bernoulli(l /2) 
and define X as 

X = B(£i, £3/3) + (1 - B){E2/2, E2/2). 

As n ^ 00, 



n 



'l-P[X^ - log ?i < x\ - log n<x^]) 



n 



-(logn+x'i)^ ( _ ^-3(logn+x2)^ _^ _ ( I _ g-2(log n+^l Ax^) 



which implies (|2.ip holds on E = [—00,00]^ \ {(—00,00)} with v (([— oo,x^] x [— oo,x^])'^) = e 
and 1^ concentrates on (—00,00] x {—00}. Thus u has zero second marginal and we seek HDA on 
E^. Note that 



nP 



X' 



n 



— - — < X ,X — > X 



.(l2|Ii+^i)\ 3(l£|Il+^2) 



_2(l2|Il+x2) 



_2(l£|Ii+xi) 



{a;2<xl} 



as n — )• 00, 



and also 



nP 



X' 



log n 



I logn 2 
> x^X'' ^ > x'' 



n 



ig-(i2|il+xl)g-3(i2|Ii+x2) ^ ig-2(l2|i+xlVx2) 



£ -2(xlVx2) 



as 71 —7- 00. 



Thus, HDA exists on E^ = [—00,00] x (—00,00] with limit measure z^'"', where z^'~'([— 00, x^] x 



(x2,00]) 



-2x' 



l|^2<^.i| and u^{{x\oc] x (x^oo]) = ie-2(^'vx ) fo^ a; € E^. In fact. 



z^'"' concentrates on the line {{x,x) : x € (—00, 00)}. 

Since i'^{{x'^,oo] x (x^, 00]) > for a; G E'~', we do not seek HDA on E°. 

3.4. Estimation methods. To estimate joint tail probabilities, we require an estimate of the limit 
measure 1^^ given in Definition 13.11 and possibly given in Definition 13.31 Estimation of i^^ follows 
the same steps as in Section [2.31 so we concentrate on estimating v^. Let, {X , Xi,i = 1, 2, • • • , n} 
be iid where the distribution of X satisfies p.ip . From ()3.ip we get [H, page 139] in M+(E'~') 



(3.11) 



1 



i=l 



i/k) Xf-f(n/k) 



e{n/ k) 



e.{n/k) 



(•) ^ v^i') {k 00, n/k 00). 



We know from (j3.2p that the distribution of X'^ characterizes {e^} and {fn} and from the iid data 
{Xf : i = 1,2, ■ ■ ■ ,n}, we can construct estimators of e{n/k) and f{n/k) denoted by e(n/k) and 
f{n/k) [H, page 93] such that 



(3.12) 



e{n/k) p f{n/k) - f{n/k) p 



e{n/k) 



e{n/k) 



0. 
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Since the limits in ()3.12p are constants, we get joint convergence in M^{K'~^) x M^, 



a-i-i^ MV". f{n/k) - f{n/k) e{n/k) \ ^ ^ ^ 

Apply the continuous mapping theorem to p.lSp using the map (z^(-),6, a) i— t- u{a[{-) + h]) to get 
in M+{¥P) 

— 1 " 

(3.14) i/n(.) : = - ^ e x^-h,./k) \ (') ^ {k ^ oo, ^ oo). 

i=l \^ e(n/fe) ' l(n/k) j 

This estimator of is non-parametric and as in Section 12. 3| we exploit the semi-parametric 
structure of by estimating 7'"' and the standardized hidden spectral measure S^. The parameter 
7^ is the extreme value index of the distribution of X'^ so estimating this from iid data {X? : 
i = 1,2,- •• ,n} is standard 0, page 65]. We obtain an estimator of S^{-) by modifying ()2.23p 
to account for the difference between E*^ and E'~'. Since the first step is to construct a consistent 
estimator of defined in (j3.6p . define 

(3.15) B"/ := I {3 : X] > X/} | and i?f := | {j : X| > } | , (i = 1, 2, • • • , n, ) 

where | • | denotes size of a set. Observe i?^'^ is just the anti-rank of Xf and thus 1 < i?^'^ < n 
for i = 1,2, ■ ■ ■ ,n. Also, < R^' < n for i = 1,2, ■ ■ ■ ,n. An estimator of is obtained from the 
convergence in M+([0, 00] x (0,oo]): 



- 1 " 

(3.16) ^ni-) ■= ]:^^{k/Rl'\ k/Rf ^ '^"(■) [k^ 00, n/k 

1=1 



00 



The verification of (I3.16P follows the steps used in the proof of Proposition l2.9l Changing coordinate 
system in (j3.16p leads to an estimator of i^i x «S'^(-) from the convergence in M+((0,oo] x [0, 00]) 

I " 

(3.17) ui X S^n{-) := - ^e(^^/^2,2^^2.2/^i,2-j(-) ^ i^i x S"{-) {k 00, n/k 00) 

i=l 

and this produces an estimator of 5^ since in M+([0,oo]), 

Ei=iefc/R2,2([l,oo]) 

This estimator may be modified as in ()2.3ip using the continuous bijection TR : [0, 00] 1— )• [0, 1] 
defined hyTR:x^ x/{l + x) to get in M+([0, 1]), 

(3.19) S^i-)oTR-'^ ^ S^'oTR-^i-). 

This summarizes how to obtain consistent estimators for extreme value index 7'^ and the stan- 
dardized hidden spectral measure 5'^. 

4. Detection of HDA 

Since HDA is a generalization of HRV, it is not surprising that the detection techniques have 
similarities to those used for HRV; see Mitra and Resnick [I4] and Resnick 15|, pages 316-340] . 
However, we deviate from the standard MEVT by assuming ()2.ip instead of ()1.2p and so proceed 
carefully. A first step is to detect the presence of asymptotic independence, which traditionally has 
been done with a density plot of a spectral measure after non-parametric transformation to Pareto 
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scale pages 316-321]. When asymptotic independence is present, we seek HDA. We consider 
how to define an appropriate spectral measure for detection of HDA. 

Define C/^^H') = 1/(1 - where iJ(^) is the distribution function of X'^^\ From I^Jij we 

get on (0, oo]^ that 

(4.1) nP [([/W(X^)/n,C/(^)(X2)/n) G •] A z>(-), 

where P(-) is a Radon measure on (0, oo]^ related to the limit measure in (12. ip by 

(4.2) P((x\oo] X (x^oo]) = 1/ ((^^(2;^),(X)] X (V'^(x2),oo]) (a; G (0, oo]^), 

and ijj is defined in ()2.4p as ^(y) := z^({z G E : y}). The measure i>(-) satisfies the scaling 

(4.3) z>(c-) = c"^P(-), c > 0, 
and convergence in (j4.ip is equivalent to 

(xW) / [/(i)(xi) ;7«(x2) 



(4.4) nP 



G • 



1^1 X 5(-), 



n ' \^C/(i)(X(i))' C/(i)(XW) 

on (0, cxo] X 5'ii.^^\ where i/i((y,oo]) = y"^ for y > 0, dii^^^ = {x € (0,oo]^ : x^^^ = 1} and 5" is a 
probability measure on The standardized spectral measure S is related to i) in (j4.ip by 

(4.5) D (^^x G (0,oo]2 : x^^^ > r,x/x^^^ G a}) = r-i5(A), r > 0, Borel set A C 6H^^\ 
To estimate this measure S{-), we define variants of the anti-ranks 

(4.6) r]'^'^ = \{j : X^'^ > Xl}\ and E^'^^ = \{j : xf^ > Xf}\ (1 < i < n). 
A consistent estimator of S{-) is obtained from the convergence in 



(4.7) 5„:= ^^..^ ^ ^ ^5(.), 

{k oo, n/k — > oo). The continuous bijection T : (JK^^^ i-^- [0, 1] given hy T : x x'^/{x^ + x"^) 
transforms g2D to Sn o T'^ ^ S o T-^{-) in M+([0, 1]). A density plot of S„ o T'^ is easier to 
analyze because [0, 1] is a nicer space than (Jb^^^^. 

Analyzing the density plot using the points of Sn o T^^ should yield evidence falling into the 
following categories: 

(i) The distribution SoT^^ concentrates near 0, so remove {{x, y) G : y = — oo} and seek HDA 
on = [— oo, oo] X (— oo, oo] or its first quadrant analogue, depending on the distribution of 
the second component of the random vector; see Remark 13.21 ([2|). 

(ii) The distribution S o concentrates near 1, so remove {{x,y) £ M? : x = — oo} and seek 
HDA on (—00,00] X [—00,00] or its first quadrant analogue depending on the distribution of 
the first component of the random vector. 

(iii) The distribution S o concentrates near and 1, so remove {(x,y) G : x = —00} U 
{(x,y) G : y = —00} and seek HDA on = (—00,00] x (—00,00] or its first quadrant 
analogue, depending on the distribution of the smallest component of the random vector; see 
Remark O®. 

(iv) The distribution S o T^^ does not have any of the above properties; we have no evidence for 
asymptotic independence and we do not consider HDA. 



HIDDEN DOMAIN OF ATTRACTION 



17 



We summarize our detection strategy and for concreteness assume S oT"^ satisfies category (i): 
Check whether X"^ belongs to some maximal domain of attraction using a Hill or Pickands plots. 
If so, conclude HDA exists on E'~'. Then consider whether HDA exists also on E*^ by examining 
a similar kernel density plot formed by using the points of the estimator of o [TR)~^ given in 

5. Conclusion 

We defined HDA as a generalization of HRV and have shown by example that for some random 
vectors, HDA exists but HRV does not. Using similar methods as in HRV, we outlined detection and 
estimation methods for HDA. These methods are given to show what is possible and to emphasize 
there is a gap that such methods can fill to provide improved estimates of probability of simultaneous 
exceedance by components of a risk vector. However, we have not implemented the methods nor 
demonstrated utility by analyzing data. This will come in the future. 

Our discussion here is restricted to two-dimensions. As observed for HRV [13], extensions to 
higher dimensions are not always straightforward and involve subtleties. In particular, in higher 
dimensions there are many more ways domains of attraction could be hidden and many more 
subspaces to explore for behavior that helps to estimate risk probabilities. 

As with HRV (l^ . our detection and estimation methods are exploratory and our estimators are 
only provably consistent. More formal statistical theory is needed to turn exploratory methods into 
confirmatory ones. 
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